Observation of a metallic superfluid in a numerical experiment 
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We report the observation, in Monte Carlo simulations, of a novel type of quantum ordered state: 
the metallic superfluid. The metallic superfluid features ohmic resistance to counter-flows of protons 
and electrons, while featuring dissipationless co-flows of electrons and protons. One of the candidates 
for a physical realization of this remarkable state of matter is hydrogen or its isotopes under high 
compression. This adds another potential candidate to the presently known quantum dissipationless 
states, namely superconductors, superfluid liquids and vapours, and supersolids. 
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At low temperatures, fluids become dominated by the 
wavelike nature of their constituent particles when the 
thermal de Broglie wavelength exceeds the interparticle 
separation. Such quantum fluids usually feature super- 
conductivity or superfluidity, which however may be de- 
stroyed by topological line-defects (vortices) threading 
the entire system. Vortices may be induced by a mag- 
netic field or rotation^ , or by thermally excited transverse 
phase fluctuations of the macroscopic wavefunction of su- 
perconductors and superfluidsSiii. In a system of rapidly 
growing interest, a two-component superconductor^, vor- 
tices yield dramatic physical consequences. Here we re- 
port the first observations, in a numerical experiment, 
of a novel type of quantum fluid originating in aggre- 
gate states of vortex matter. Increasing temperature, a 
composite vortex lattice melts into a composite vortex 
liquid, whence superconductivity is lost while superflu- 
idity is retained, yielding a metallic (ohmic) superfluid. 
At higher temperature, another unusual transition oc- 
curs where the composite vortex liquid "ionizes" into a 
"plasma" of constituent vortices, destroying superfluid- 
ity. 

Recently, there has been considerable interest in the- 
ories of superconductors with several superconducting 
components coupled by a magnetic field, but with no 
possibility of Josephson tunneling of one component 
into another. This is predicted to occur in a wide va- 
riety of physical systems, most notably in condensed 
matter such as hydrogenic atoms subjected to extreme 
pressure^-^-^-^-^-^^ or effective theories of easy-plane quan- 
tum antiferromagnetsii. Renewed interest in the long 
sought liquid metallic hydrogen (LMH) is due to recent 
ab initio calculations^^ along with a breakthrough in syn- 
thesis of ultrahard artificial diamonds, essential for ob- 
taining the required extreme pressures in anvil cells 
These facts, along with a recent measurement of an un- 
usually low melting temperature of dense Na^^, strongly 
hint at a realization of LMH in the not too distant future. 
Thus, understanding the superconducting properties of 
this system is important, since magnetic field experi- 
ments can be conducted in high pressure anvil cells, possi- 
bly confirming the realization of this novel state of matter 
in a terrestrial laboratory. It is commonly accepted that 



LMH is abundant in the interior of Jupiter and Saturn 
and quite possibly also present in some of the known 200 
extrasolar giant planets^^. In these cases, however, LMH 
is conjectured to exist in the classical metallic regime at 
several thousand degrees Kelvin. In contrast, the state 
projected to exist in Refs. 6 is a ground state quantum 
fluid, and is the one we focus on in this Letter. 

Realization of LMH could well constitute the next 
milestone in quantum fluids. It is projected to fea- 
ture Cooper-pairs of both electrons and protons at low 
temperatures^. The resulting quantum fluid differs rad- 
ically from previously known quantum fluids, in that its 
physical properties cannot be classified exclusively as a 
superconductor or a superfluid^. Remarkably, such a 
system features both superconductivity and superfluidity 
which appear as collective phenomena corresponding to 
CO- and counter-flows of two species of Cooper-pairs, with 
a complicated interplay between them. Thus, LMH has 
been conjectured to sustain phase transitions connect- 
ing a superconducting and superfluid state to a metallic 
state featuring a superfluid mode^, or to a superconduct- 
ing state with no superfluid mode. 

The transition from a state featuring superconductiv- 
ity and superfluidity to a state where superfluidity is lost 
and superconductivity is retained, has recently been ob- 
served in a large scale Monte Carlo (MC) simulationiS. 
However, the remarkable possibility of a transition from a 
"composite" vortex lattice (superconducting state where 
vortex matter forms a solid) to a "composite vortex liq- 
uid" (metallic superfluid state of the system) along with 
a subsequent transition from a composite vortex liquid 
to "vortex line plasma" , has thus far not been confirmed. 
It is the purpose of this Letter to report an observation 
of the metallic superfluid, not dealt with in simulations 
previously, in a numerical experiment. 

The superconducting phase of LMH is given by the 
Ginzburg-Landau model with two scalar fields ^Q^^(r) 
and ^Q^^(r) describing superconducting condensates of 
protons and electrons, respectively. It is defined via the 
energy density 
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Here, M^'^^ and M^^^ are the masses of the condensates, 
the covariant derivative is given by D = V — ieA(r), 

and V({|^Q"^(r)|}) is the potential term. The par- 
ticular form of V({|^Q"^(r)|}) is not essential for the 
large-scale physics we address, but its dependence only 
on l^o^^l reflects the fact that Cooper-pairs of elec- 
trons cannot be converted into Cooper-pairs of protons, 
and vice versa. For the issues discussed in this pa- 
per, it suffices to work in the phase only approximation 
^[)''^(r) = |^[)"^|exp[il9^")(r)] where |^[)"^| is treated as 
a constant In LMH, the two order parameters cor- 
respond to electronic and protonic Cooper-pairs. More- 
over, Eq. (pP) may be rewritten as follows^iSiiS 
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The neutral and charged modes, described by the second 
and first term in Eq. are explicitly identified. 

The topological objects of Eq. (QJ are vortices of type- 
1 and type-2 defined by a 27r winding in O'^^^ and 6>'^^\ 
respectively. The interaction potential between these 
vortices is a superposition of a Coulomb potential and 
a Yukawa potential, arising out of the neutral and the 
charged mode, respectively. (For a derivation of the for- 
mulae for the vortex interaction, see Ref. 10). The en- 
ergy of a vortex associated with ±27r winding in 0^^"^ —0^'^^ 
is logarithmically divergent. As the neutral mode tends 
to lock V^^^^ to V6^'^\ in order to minimize the second 
term in Eq. the vortices of type-1 and type-2 pair 
up into a composite vortex for which V(^^^^ — 0^'^^) = 0. 
Therefore, a composite vortex is an object where a type- 
1 and a type-2 vortex are co-centered and co-directed in 
space, and the Coulomb part of the pair-potential exactly 
cancels. The screened potentials add, but the associated 
overall energy is finite^'^'^'^^. In the presence of an ex- 
ternally applied magnetic field, the ground state of the 
system is a lattice of co-centered vortices of type-1 and 
type-2, a composite vortex lattice. 

We have performed MC simulations on Eq. (^) at finite 
temperature using local Metropolis updating on the fields 
^^""^^(r), ^^^^(r), and ^(r). The system size we have used 
is L X L X L, with L = 120. The coupling constants inves- 
tigated are IV^^^^P = 0.5, IV^^^^p = 1.0 and e = 1.0 and 
the external magnetic field isB = VxA(r) = (0,0,27r/) 
with / = 1/20. Thus, there are 20 plaquettes in the xy- 
plane for each field-induced vortex. The external mag- 
netic field is imposed by splitting A(r) = Ai?(r) + Ao(r), 
where A^? is free to fluctuate subject to periodic bound- 
ary conditions, and Aq = (0, 27rx/, 0) is kept fixed. We 
have chosen the amplitude ratios |?/^^^^ p = 2.0 for 

numerical convenience. We emphasize that the results 
in this paper are dictated by symmetry, and the physi- 
cal picture we present will thus be representative also for 
LMH, where |2/|^(i) |2 _ ^^s 



To find the lattice ordering of vortices we compute the 
planar structure function 6'^"^(k^) of the local vorticity 
n(")(r), defined by A x [Al9(") - eA] = 27rn(")(r), given 

by5M(k^) = (lEr e^^— p)/(/L3)2.Here, 

is the lattice difference operator, A^O^^"^ — eA^ G [0, 27r), 
r runs over the possible positions of the vortices, and 
and are perpendicular to B. If vortices form a lattice, 
^(«)(k^) win exhibit a six-fold symmetric Bragg struc- 
ture, and feature a ring-structure in the vortex liquid 
phase. 
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FIG. 1: (Color online) Results of numerical experiments on 
a two-component vortex system, for IV^*^"^-*!^ = 0.5, IV^*^^-*!^ = 
1.0, and e = 1.0. Panel a: Structure functions for protonic 
S^'^\'K) (red) and electronic S**^^-* (K)(blue) vortices, for K = 
(7r/4, 27r/5). They drop to zero discontinuously at Tm where 
the system looses its superconducting properties, but retains 
superfluidity, as evidenced by a finite helicity modulus Tx 
(black line). Tx serves as an order parameter in the vortex 
liquid phase, dropping to zero at Tc. Panel b: A^co (green) is 
finite across the melting transition, but drops continuously to 
a small value at Tc, where it has a kink. At Tc, the specific 
heat Cv has an anomaly. 

The MC results are given in Fig. Q showing the 
structure function 6'^^^(K) (red) for protonic vortices, 
and 6''^^^(K) (blue) for electronic vortices, where K = 
(7r/4, 27r/5) is a Bragg vector. In the low temperature 
regime both functions are finite, but decrease gradually 
as the temperature is increased. At Tm, 6'^^^(K) and 
6''^^^(K) vanish at the same temperature even though the 
ratio of the bare stiffnesses of the condensates is 2.0. 
Moreover, both structure functions vanish discontinu- 
ously, the hallmark of a first order melting transition of 
the composite vortex lattice. 

To probe the splitting of the co-centered vortices 
into constituent vortices, we compute the vortex co- 
centricity, defined as A^co = ^co"^ — ^c^^ where Net'' = 
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if i = j and zero otherwise. Therefore, A^co is the frac- 
tion of type-2 vortex segments that are co-centered and 
co-hnear with type-1 segments. We find that passing 
through the first-order melting transition renders Nco un- 
affected, whence we conclude that the observed transition 
is a melting of a composite vortex lattice. 

In one-component type-II superconductors, melting of 
the vortex lattice amounts to a complete destruction of 
dissipationless currents'^. To follow the fate of the super- 
fluid mode of the two-component system in the vortex 
liquid state, we measure the ordering in 7(r) = ^^-'^^(r) — 
^^^^(r). To probe the global phase coherence in this vari- 
able we consider the helicity modulus T^, equivalently 
the superfluid density, given by = [(c) — 77 (s^)] /L^, 
where c = J^lr ^os (A^7(r)), s = sin (A^7(r)), and 
rj = /3|?/;(i)p|^(2)|2/2^2 When the composite vortex 
lattice melts, destroying superconductivity^, the super- 
fluid density remains unaffected^ cf. Fig. ^ Hence, this 
transition separates a superconducting superfluid from a 
metallic superfluid state. 

Thus, we have identified a transition from a composite 
vortex lattice into a composite vortex liquid. Increas- 
ing the temperature further, we find a phase transition 
between the composite vortex liquid and the 'Honized vor- 
tices'^ plasma^ as evidenced by the vanishing of the su- 
perfluid density T^^, cf. Fig. C] We observe that this 



transition is accompanied by a pronounced anomaly in 
the specific heat Cy , indicating that this is a critical phe- 
nomenon. This is corroborated by the following physical 
argument. The liquid state of vortices in an ordinary su- 
perconductor is a state where translational symmetry is 
restored and superconductivity is lost^. In the composite 
vortex liquid state, every electronic vortex is accompa- 
nied by a protonic vortex performing only finite excur- 
sions away from the electronic vortex line. Therefore, for 
every plane slicing an electronic vortex in a direction per- 
pendicular to the magnetic field, it is possible to identify 
a finite length closed contour which also encompasses an 
accompanying co-directed protonic vortex. Along such 
a contour, there is no nontrivial winding in the phase- 
difference 7. Therefore, melting of the composite vortex 
lattice into a composite vortex liquid does not restore the 
broken global /7(l)-symmetry associated with 7. On the 
other hand, in the vortex plasma state one cannot find 
a protonic vortex accompanying every electronic vortex, 
implying a disordering of 7. Hence, the associated global 
U (l)-symmetry is restored during the "vortex ionization" 
transition taking place within the vortex liquid, and it is 
therefore in the 3Dx?/ universality class. 

To gain further insight, we have also extracted vor- 
tices and visualized snapshots of configurations of vortex 
matter in small segments of the system. The results are 
shown in Fig. [21 



Panel a of Fig. [21 shows the side- and top- views of vor- 
tex matter at a small, but finite, temperature when two 
species of vortices perform only small excursion from each 
other. The insets show the structure function in momen- 
tum space of the protonic (thin red) and electronic (thick 
blue) vortices. Parallel and co-directed vortices inter- 
act with each other as positively and negatively charged 
strings, and the splitting is a temperature-induced fluc- 
tuation. Panel b shows snapshots of the vortex matter 
when the system is heated above the temperature of vor- 
tex lattice melting. The protonic (red) and electronic 
(blue) structure functions have developed ring-like struc- 
tures, characteristic of a vortex-liquid in both the pro- 
tonic and electronic sectors. While the vortices perform 
stronger excursions from each other, these excursions are 
still limited and one can always identify a red line at- 
tached to any given blue line. Thus, co-centricity of 
protonic (red) and electronic (blue) vortex lines is still 
largely intact. In such a configuration, a dissipation- 
less electrical current cannot propagate in any direction. 
Quite remarkably, however, co-directed currents of pro- 
tonic and electronic pairs can propagate through this sys- 
tem without dissipation, as evidenced by the measure- 
ments of the helicity modulus given in Fig. ^ Panel 
c shows a state of vortex matter which occurs above 



the vortex "ionization" temperature. The co-centricity 
of vortices is strongly reduced, see Fig. 

This is also reflected in the subtle difference between 
protonic (red) and electronic (blue) structure functions 
in passing from panel b to panel c. While the struc- 
ture function of the electronic vortices essentially is un- 
affected by passing through the temperature Tc, the struc- 
ture function for the protonic vortices is distinctly further 
isotropized. Namely, the relative increase of the uniform 
background for the structure function of protonic vortices 
is a manifestation of the fact that protonic vortices suffer 
a vortex-loop proliferation transition inside the metal- 
lic vortex liquid phase. We have argued above that this 
transition belongs to the 3Dxy universality class, i.e. the 
same universality class as the superfluid-to-normal fluid 
transition in liquid helium ^He^^. 

Concluding, we report the observation, in a numeri- 
cal experiment, of a novel dissipationless quantum state 
of matter, namely the metallic superfluid. Such a state 
might be realized in hydrogen or deuterium, if those sys- 
tems were to take up a projected low-temperature liq- 
uid metallic state at an extreme pressure. The high- 
est pressure obtained to date appears to be around 
320GPa^^. However, recent breakthroughs in artificial 
ultrahard diamond synthesis technology^^ represent sig- 
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FIG. 2: (Color online) Snapshots of the states of vortex matter and momentum space structure functions generated from MC 
simulations, taken at three different temperatures: T = 0.50 (T < Tm), T = 0.72 (Tm < T < Tc), and T = 0.86 (T > Tc). The 
snapshots are extracted from a small segment (15 x 15 x 15) of the vortex system. Each thick frame contains a sideview and a 
topview of the vortices, as well as the protonic structure function 5'*^^'*(k^) (red) and the electronic structure function S^'^\'k±) 
(blue). Panel a: For T < Tm the vortices are arranged in a co-centered lattice. Electronic (thick blue) and protonic(thin red) 
vortices only perform small excursions from each other. Both structure functions exhibit a sixfold Bragg pattern characteristic 
of a vortex lattice. Panel b: For Tm < T < Tc the composite vortex lattice has melted as illustrated by the ring pattern 
in the structure functions. The electronic and protonic vortices perform stronger excursions from each other, but essentially 
remain co-centered. This is the superfluid metallic phase in which co-directed currents of protonic and electronic Cooper-pairs 
can propagate without dissipation. Panel c: For T > Tc the superfluidity is lost and the electronic and protonic vortices are 
no longer co-centered. The proliferation of the protonic vortices is reflected in the increase in the uniform background of the 
protonic structure function. 



nificant progress towards achieving extreme pressures in 
diamond anvil cells. Thus, a metallic superfluid might 
be the next "super" state of matter to be realized in the 
laboratory. 
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